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Abstract 



We formulate the calculation of the ground-state wavefunction and energy 
of a system of strongly correlated electrons in terms of scattering matrices. 
A hierarchy of approximations is introduced which results in an incremental 
expansion of the energy. The present approach generalizes previous work 
designed for weakly correlated electronic systems. 
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I. INTRODUCTION 



A microscopic wavefunction-based description of electron correlations in 
the ground state of extended systems, i.e., large molecules or solids remains 
a challenging problem. This holds particularly true when the correlations are 
strong. Much progress has been achieved though over the last twenty years. 
For example, it has become clear that local operators have to be applied for de- 
scribing the correlation hole of the electrons [1-6] . Otherwise calculations with 
controlled approximations for extended systems become nonfeasible. When 
starting from a self-consistent field (SCF) wavefunction for the ground state 
denoted by | <3?o > the wave operator Q, which transforms it into the true 
ground-state wavefunction | ty$ > must be constructed from local operators 
A v . They describe one-, two- or generally multiparticle excitations out of 
| $o >• The fact that the A v are local operators implies that the creation 
and annihilation operators appearing in A v refer to local orbitals instead of 
canonical orbitals or Bloch states. 

The first ground-state calculations based on the use of local operators were 
done for diamond [7]. However, it was not until the method of increments was 
pointed out [5,8], that a large number of solids were successfully treated by 
quantum-chemical configuration- interaction (CI) techniques (see, e.g., [8-10]). 
Although most of the calculations were done with respect to the ground-state 
wavefunction it was shown that the same ideas can be also applied to the 
calculation of excited states, i.e., energy bands of a solid [11]. 

The calculations described above were done for systems in which electron 
correlations are not too strong, i.e., for which the SCF configuration is a good 
starting-point. If the electronic correlations are strong one should not start 
from the independent-electron approximation when attempting to calculate 
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the true ground-state wavefunction. It has been suggested that the method of 
increments can be also applied to strongly correlated systems by performing 
multi-configuration SCF (MCSCF) or complete-active-space SCF (CASSCF) 
calculations for localized orbital groups [12]. But the theoretical formulation 
of such an approximation scheme as well as the form of the ground-state wave 
function have remained unclear. 

The aim of the present communication is to provide a basis for the incre- 
mental method when the correlations are strong. For that purpose we have to 
modify a derivation of an incremental scheme given in Ref. [13]. We want to 
show how one can construct a cumulant wave operator for strongly correlated 
electron systems. This operator defines the ground-state of the system. 

The paper is organized as follows. In the next section a frame is provided 
for the computation of the ground state of a strongly correlated extended 
electron system by means of quantum-chemical methods. Section III demon- 
strates explicitly how the calculations have to be done in practice. Section IV 
contains a summary and the conclusions. 



II. FORMALISM 



Starting point is the Hamiltonian H of the electronic system 

ija ijktao 1 

It refers to a given basis set fi(r) for which usually Gauss- type (GTO) or 
Slater-type orbitals are chosen. The matrix elements Uj and Viju refer to the 
one-electron (including kinetic energy) and two-electron interaction energy, 
respectively. We split H into a self-consistent field part Hscf and a residual 
interaction part H res . The latter is given, for the closed-shell case, e.g., by 
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H, 



res — 



ijki L era' 




The one-particle density matrix Pij is defined by 



< a t* a 3° > 



a 



and < ... >=< &SCF I ••• I &SCF > with | &SCF > denoting the ground 
state of Hscf- It is noticed that H res consists of a constant term plus one- 
and two-particle excitations. We subdivide H res further into 



Here Hi denotes that part of H res where excitations are restricted to a 
group of localized orbitals belonging to center / (e.g., an atom or a bond), 
while Hu contains the part in which the centers / and J are involved in 
the process of creating holes. For example, van der Waals type interactions 
between atoms (or bonds) / and J belong to Hjj since they involve a one- 
particle excitation on each of the atoms (bonds) / and J. This division of H res 
differs from the one in Ref. [13] where H res was decomposed into contributions 
corresponding to different pairs of holes generated out of the self-consistent 
ground state | &scf >■ 

We characterize the ground state of H by the cumulant wave operator 
I f2) [14, 6]. The conventional wave operator Q is defined as the one which 
transforms | <&scf > into the exact ground state | ^0 >, i- e -> I ^0 >= ^ I 
&SCF >• The cumulant wave operator | ft) differs from Q, in that it is defined 
in a space with a different metric than that, e.g., of the conventional Hilbert 
or operator space. Here the metric is defined by the following bilinear form 
of two general operators A and B. 



H res = Y Hl + Y HlJ - 



1 1, J 



(A I B) = < $scf I A+B I $scf > c • 



The upper script c implies taking the cumulant of that expression [14,15,6]. 
The usefulness of cumulants in quantum mechanics was stressed by Kubo 
[15] who generalized earlier work of Ursell and Mayer in classical statistical 
mechanics [16]. Cumulants have the advantage that they correspond to linked 
clusters. Therefore the problem of size consistency (or extensivity) does not 
exist when quantities are expressed in terms of them. For further information 
we refer to, e.g., Ref.[6]. In accordance with the above the ground-state energy 
is given by 

E = (H\Q). 
The cumulant wave operator is of the form [6] 



0,) = lim 



z — H 



In accordance with Ref. [13] we define a scattering operator 

In proceeding we use an argument developed in [13] which can be consid- 
ered as a generalization of Faddeev's equation. Faddeev derived an equation 
[17] which expresses the scattering operator of a three-particle system in terms 
of the scattering matrices for the different two-particle channels. For the latter 
often analytic solutions can be found. More generally, we aim at expressing 
the scattering operator of an N electron system in terms of scattering matrices 
of simpler subsystems. With this goal in mind we introduce the operators 

A n = lim \- Hi , A u = lim \- H u . 

z^O z - H S CF 2 ->o z - H S CF 



With their help we rewrite 

' N 



s) = E 



n=l 



where N is the number of centers, i.e., atoms (bonds) in the system. We 
introduce Greek letters to denote pair labels I J so that Yl, I M.J >= J2 \ A a > 

IJ a 

and decompose | S) into terms with one Greek index and a remaining part. 
Therefore we write 

* = E (e^ + E^ 

a \n=l J 

= E s <* + E 

The T Q/ g are defined as follows: When (10) is decomposed into different 
terms we include in T a p all those which begin with A a from the left followed 
by Ap as the first factor different from A a . For example, terms of the form 
A a Ap..., A a A a Ap..., A a ApAp... are all included in T a p. With this definition 
we can write 

T a p = (A a Aj3 + A a ApA a + ...) • 

• (i+ E ^+ E T ~tX 

V 1^<*,P;5 ) 

When adding T a p + Tp a one notices that the first bracket is nothing else 
but the scattering operator S a p of a Hamiltonian H$cf + H a + Hp, except 
that the contributions S a + Sp are missing. Therefore we can write 

T a p + Tp a = (S a p - S a - Sp) 11+ y~] S 7 + ^2 T-yS J ■ 



or 

S = E 5 ° + E (SaP - S a - Sp) ll + £ 5 7 + X! T 7<5]- 

Here < a (3 > denotes different pairs. 

One expects that the matrix elements of Hj are generally much larger than 
those of Hjj. Therefore it seems advantageous to introduce approximations 
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to I S) by resumming the right hand side of the last equation according to 
different numbers of sites involved. To lowest order | S) is therefore given by 

\S) = J2Sj, (14) 
i 

where Sj (= S a with a = II) is the scattering operator of Hamiltonian 

H SC F + Hi- 

In this single-site approximation 

E = E SCF + (H\S) (15) 

= E S cf + ^2 e i 
I 

with ej = (H | Si). In next order we include the terms S a with a = IJ 
and T a p with a and (5 being iY, J J, IJ and .11. In this two-sites approximation 
the second bracket in (12) is replaced by unity. By adding up the different 
contributions we find 

I S) = I s i) + E I S " ~ S i ~ S J) 

i <u> 

where Sjj is the scattering operator belonging to Hscf+Hi+Hj+Hij+ 
Hji- 

Furthermore, 

E = E SC F + E e i + E e u ( 16 ) 

/ <u> 

with ejj = (H | Su) — e/ — ej. This procedure can be continued. By 
including also three-sites terms we find 

I s) = ]T | Sj) + £ | ss u ) + £ I as/./*) + ... (17) 

/ </J> <I.JK> 

with | (JS'/j) =| Sjj)— I Si)— I 5"j) as before and 
I SS IJK ) =| I SSu)- | <J5ja:)- | SS ki )- \ Sj)- \ Sj)- \ S K ). (18) 
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Here | Sijk) is the scattering operator of a Hamiltonian H = Hscf + 
Hi + Hj + H K + Hu + H IK + H JK + Hjj + H KI + H KJ . This approximation 
scheme can be continued until finally the exact scattering operator | S) is 
obtained. The corresponding energy expression (15) is then 

Eq = E S CF + Y €l + Y ElJ+ Y e UK + -~ 
I <IJ> <IJK> 

with 

eijK = (H | Sijk) — eu - cik - £jk — e/ — ej - ex 

etc. 

The advantage of the above formalism is that we have reduced the ground- 
state calculations for extended systems like solids to the computation of single- 
center, two-center etc. scattering matrices. These matrices can be determined 
by means of quantum-chemical methods whereby all the other electrons in 
I ®SCF > are kept frozen. There is no difficulty in calculating the scatter- 
ing matrices also when the electrons are strongly correlated. For example, 
multiconfiguration SCF (MCSCF) calculations or complete-active-space SCF 
(CASSCF) calculations, followed by a multireference CI (MRCI) treatment, 
serve that purpose. Strong correlations can therefore be treated with a high 
degree of accuracy in ground-state calculations for solids. 

III. APPLICATIONS 

In the following we want to outline in some more detail how ground-state 
calculations for strongly correlated electron systems have to be performed. We 
limit ourselves to insulators or semiconductors, i.e., systems with a gap in the 
excitation spectrum. Starting point is the Hamiltonian (1) acting in a space 



spanned by a properly chosen basis set of GTO's. After a SCF calculation 
has been performed, e.g., by using the program package CRYSTAL [18] or 
the code developed by Shukla et al. [19], one has to express the SCF orbitals 
in the form of orthogonal localized Wannier orbitals. This is achieved either 
by an a-posteriori localization procedure if CRYSTAL is used, or by using 
Shukla's program which yields directly the occupied SCF orbitals in localized 
form. The corresponding creation and annihilation operators are denoted 
by cf a ,C£ a . For the strongly correlated electrons in the system, e.g., the 
d electrons of a transition metal atom or the / electrons of a rare earth or 
actinide atom, we need to express in localized form not only the occupied, but 
also the unoccupied (virtual) d or / orbitals. Finding them poses no problem 
and we include them in the set of operators cf a , C( a . In a next step the residual 
interactions H res are expressed in terms of the c^.,q ct operators. Except for 
the strongly correlated electrons we express only the annihilation operators in 
(2) in terms of the q ct . (The external orbitals have to be localized, but they 
need not being orthogonal to each other, cf. e.g. [20].) The strong correlations 
are treated by a CASSCF calculation. The choice of the active space depends 
on the scattering matrix we want to calculate. For the determination of the Si 
we choose the strongly correlated orbitals of that center (e.g., the localized d 
or / orbitals of an atom or the bonding and anti-bonding orbitals of a bond) 
for the active space. With all electrons kept frozen except those on center 
I, the CASSCF calculation with the active space involving orbitals on that 
center only accounts for the strong intra-atomic (or intra-bond) correlations. 
When Si j is calculated we must include in the active space the localized d(f) 
orbitals (or bonding/antibonding orbitals) on centers I and J. The resulting 
ground state is denoted by | >• After every CASSCF calculation, the 
remaining (weak) correlations may be accurately taken into account by means 
of calculations within some variant of the coupled electron pair approximation 
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(e.g., MRCEPA or MRACPF with single and double excitations into the 
external space). In case that the two-center scattering matrices Su are not 
sufficient for an accurate determination of | f2) one can extend the calculations 
to three-center scattering matrices Sjjk- 

Let us apply now the above considerations to simple molecular exam- 
ples involving hydrocarbons. We address the question of how to properly 
treat such systems in the strong-correlation limit of dissociation into separate 
atoms; the ultimate goal of our investigation being a unified treatment of Coo 
(diamond) for a wide range of C-C internuclear distances. Specifically, we per- 
formed calculations for CH4, C2H6, C3H8, and C5H12 (neopentane); all angles 
were kept fixed at 109.47°; standard equilibrium C-H and C-C bond lengths 
of 1.102 and 1.544 A , respectively, were uniformly scaled by factors /, with 
1 < / < 100. We concentrated on strong correlations, i.e., restricted our basis 
set to the single-zeta level ((9s4p)/[2slp] for C, (4s)/ [Is] for H, using subsets 
of Dunning's correlation consistent valence double-zeta sets [21]). Since SCF 
calculations, even with such small basis sets, encounter severe convergence 
difficulties for large internuclear distances (and yield physically unreasonable 
energies anyway, high above the dissociation limit), we started from local- 
ized two-center orbitals generated as follows. We combined sp 3 hybrids on 
the C centers with each other and with H Is orbitals to form bonding and 
anti-bonding LMOs (coefficients ±1) between next-neighbour atoms; we then 
Gram-Schmidt-orthogonalized all valence orbitals to the C Is cores, symmet- 
rically orthogonalized the bonding LMOs among each other, proceeded by 
Gram-Schmidt-orthogonalizing the anti-bonding LMOs to the bonding ones, 
and finally symmetrically orthogonalized within the anti-bonding space. With 
this construction, we can build up an SCF-like closed-shell state, with all 
bonding LMOs doubly occupied, which should resemble the true SCF ground 
state of Coo, in the limit of large cluster size (and in a minimal-basis set rep- 
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resentation, of course). We then defined groups of orbitals, pairing each of 
the bonding LMOs with the corresponding anti-bonding one, and performed 
CAS CI calculations [22 - 24] with one of the groups active in turn — this 
leads to correlation-energy increments Aecc and Aecn describing the break- 
ing of a CC or a CH bond in a frozen closed-shell environment. Their sum 
provides us with a first approximation to the correlation energy of the system, 
but still not a very good one, since reorganization at the C atoms (leading 
to 3 P ground states in the limit of R — > oo) is not taken into account. We 
therefore introduced, in the next step, an atomic correction by correlating 
simultaneously all eight LMOs (bonding and anti-bonding ones) related to 
a given C atom, in a CASCI calculation. This defines an atomic increment 
AAec = ec — J2i Aecx, where ec is the correlation energy of the calculation 
just mentioned, and the Aecx are the single-bond increments of the neigh- 
bouring atoms X = C, H. Again, by adding up all AAec contributions, we 
obtain an improved estimate for the energy of the system. (The next step - 
which we did not perform any more - would be to determine non-additivity 
corrections for pairs, triples etc. of atomic increments. In the limit of n-tuple 
corrections (n — > oo), this should lead to the full-CI energy of the system, 
irrespective of the starting-point chosen, i.e., irrespective of the fact that we 
did not start from a variational SCF wavefunction.) 

In Table 1, the so-obtained correlation-energy estimates for CH 4 and C2H 6 
are compared to reference values (full valence CI results). The errors of the 
'SCF' energies are huge, of course (~ 1.5 a.u. for / = 100 and ~ 0.13 a.u. 
for / = 1, in the case of CH4); about half of the error for / = 1 is due to 
the non-self-consistent preparation of the ground state. Including the bond 
correlations, Aecc and Aecff, errors are reduced by nearly an order of mag- 
nitude (to ~ 0.25 a.u. for / = 100 and ~ 0.02 a.u. for / = 1, in the case 
of CH4). The atomic correction, AAec, corresponds to full CI, in this case, 
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so we need not discuss it further. It is interesting to note, however, that it 
also yields the 'exact' result for C2H6 at large distances (i.e., the atoms are 
properly decoupled to separate ground-state entities), and it deviates by only 
1 mH from the full CI value for C2H6 at / = 1 (using the same basis set). 
Note that the error of a standard CCSD(T) calculation is of the same order 
of magnitude, at that internuclear distance, and substantially increases for 
/ > 1 (~5 mH for /=1.5). 

Table 2 gives a compilation of bond increments, Aecc ; an d atom incre- 
ments, AAec, for various hydrocarbon molecules, again for a large range of 
internuclear distances. It is seen that the Aecc are fairly stable in various 
environments; although their absolute value changes by more than two orders 
of magnitude, the maximum relative change is 3% between C2H6 and C5H12. 
The AAec are less transferable: they are invariant, of course, for / = 100 as 
they should, but for / = 1 the change from four H neighbours (in CH4) to a 
purely C-atom neighbourhood (for the central atom in neopentane) enhances 
the AAec by nearly a factor of 2. Assuming that changes in the second- 
nearest neighbour shell do not appreciably modify AAec any more, we can 
make an estimate for the infinite solid, on the basis of our results. We predict 
the correlation energy of diamond, Coo, for our single-zeta basis set, to be 
~ 4Ae C c(C 5 H 12 ) + 2AAec(C 5 #i 2 ) per unit cell, and obtain -1.8049, -.9921, 
-.4681, and -.1286 a.u., for / = 100,2, 1.5, and 1, respectively. 

It is clear that dynamical correlation effects left out in our example, have 
significant influence on the properties of diamond, as shown in our previous 
work [5,8]. Therefore, we plan to include such effects, in the future, at the 
MRCI (or rather MRACPF) level, into our calculations. This would enable a 
reliable description of the diamond potential-energy surface up to quite large 
internuclear separations. 
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IV. SUMMARY AND CONCLUSIONS 



We have shown that with the help of multicenter scattering matrices the 
cumulant wave operator | O) can be constructed even when the electron cor- 
relations are strong. The operator Q defines the exact ground state. It follows 
that the corresponding ground-state energy can be calculated in form of incre- 
ments as previously done for weakly correlated electron systems. The applica- 
tions described in Section III assumed insulators or semiconductors, because 
in that case orthonormal localized SCF orbitals can be easily constructed. In 
principle, however, the scattering matrix approach can be also formulated for 
nonorthogonal local orbitals. Within the theoretical framework outlined here 
accurate ground-state wavefunction and energy calculations become feasible. 
It remains a challenging problem to extend the theory to excited states. 
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TABLES 

TABLE I. Energies of an SCF-like initial wavefunction without/with subsequent correlation 
corrections, Ae/j and AAe^- (cf. text for definition), as a function of the bond- length scaling 
factor /, in comparison to full-CI calculations (FCI). All energies in Hartree. 



a) CH 4 



/ 


SCF 


+ EAe/j 


+ EAAe K 


FCI 


100 


-38.175213 


-39.446762 


-39.696730 




2 


-39.260582 


-39.600272 


-39.760294 




1.5 


-39.727852 


-39.894261 


-39.974412 




1 


-39.990677 


-40.100800 


-40.122505 




b) C 2 H 6 


/ 


SCF 


+ EAe,j 


+ £AAe*r 


FCI 


100 


-75.661403 


-77.894968 


-78.394909 


-78.394909 


2 


-77.493718 


-78.139308 


-78.485026 


-78.486933 


1.5 


-78.345506 


-78.647679 


-78.833566 


-78.833278 


1 


-78.882457 


-79.059017 


-79.120651 


-79.121427 
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TABLE II. Correlation corrections for C-C bonds, Aecc, an d atomic re-coupling, AAec, for 
various hydrocarbon molecules, as a function of the bond-length scaling factor /, cf. text. All 
energies in Hartree. 



a) Ae C c 


f 


C2H6 


C3H8 


C5H12 


100 


-.326241 




-.326241 


-.326241 


2 


-.138377 




-.138113 


-.137592 


1.5 


-.054757 




-.054566 


-.054178 


1 


-.012927 




-.013081 


-.013386 


b) AAe c 


/ 


CH 4 


C 2 H 6 


C 5 H 12 


100 


-.249968 




-.249970 


-.249970 


2 


-.160022 




-.172859 


-.220869 


1.5 


-.080151 




-.092944 


-.125696 


1 


-.021705 




-.030817 


-.037532 
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